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5.1 INTRODUCTION 

In ocean*color remoie sensing, approximately 9D % of the flux at the senior originate?; from atmospheric 
scattering, with the water* leaving radiance contributing the remaining 10 % of the total flux Consequently, errors 
in the rr>eaxured top*of*the atmosphere radiance are magnified a factor of 10 in the determination of water*leaving 
radiance. Proper characterization of the atmosphere is thus a critical part of the analysis of ocean*color remote 
sensing data. It has always been necessary to calibrate the ocean<olor satellite sensor vicariously, using in situ, 
ground'based results, independent of the status of the pre-flight radiometric calibration or the utility of on-board 
calibration strategies (Gordon 1998). Because the atmosphere contributes significantly to the measured flux at the 
instrument sensor, both the instrument and the atmospheric correction algorithm are simultaneously calibrated 
vicariously. 

The Marine Optical Buoy (MOBY) (Clark et al. 1997; Clark et al. 2002a; Clark ei al. 2002b), deployed in 
support of the Earth Observing System (EOS) since 1996, serves as the primary calibration station for a variety of 
ocean<olor satellite instruments, including the Sea*viewing Wide Field*of*view Sensor (SeaWiFS) (Barnes et al. 
2000). the Moderate Elesolution Imaging Spectroradiometer (MODIS) (Esaias et al. 1998), the Japanese Ocean 
Color Temperature Scanner (OCTS) (Isaacman et al. 1999), and the French Polarization and Directionality of the 
Earth’s Reflectances (POLDER) (Deschamps ei at. 1994). MOBY is located off the coast of Lanai, Hawaii The 
site was selected to simplify the application of the atmospheric correction algorithms (Clark et al. 1 997). Vicarious 
calibration using MOBY data allows for a thorough comparison and merger of ocean-color data ftom these multiple 
sensors (Wang et al. 2002) 

MOBY uses an instrument known as the Marine Optical System (MOS) to detect radiation over the spectral 
range from 350 nm to 955 nm. The MOS system contains two single*grating spectrogr^hs, a blue spectrograph 
(8SG) to measure light in the near ultraviolet and visible from 340 nm to 640 nm and a red spectrograph (RSG) to 
rrwasure light in the red and near infrared from 550 nm to 955 nm (Clark ci al. 2002a; Clark et al. 2002b). MOS 
resides in the MOBY instrument bay located at the bottom of the buoy. It is connected by optical fibers to radiance 
and irradiance pons on the three MOBY arms (denoted Top, Mid, and Bot), located at different ocean depths 
(typically l.S m, 5 m and 9 m), as well as to a surface irradiance port. MOBY measures upwelling radiance, L^, as 
well as the down* welling irradiance, E^. A fiber^iptic multiplexer in the MOBY instrument bay selects which 
spectrum is acquired by MOS. As described in Clark ei al. (2002b), ihe?;e data are used to determine the water* 
leaving radiance, 

MOBY buoys are typically deployed for 3 months lo 4 months, then retneved for servicing and repair. 
Deployments are numbered sequentially. One MOS instrument, MOS204, is used for even*buoy deployments while 
a separate instrument, MOS205, is used for odd deployments The buoys are calibrated before and after deployment 
using sources traceable lo radiorr>etTic sbndards maintained al the National Institute of Standards and Technology 
(NIST) (Clark er al. 2002a) During a deployment, the buoy's radiometric stability is checked using on*board 
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calibraiion M^urces and monthly lamp calibrations by divers (Clark ct al 2002a). The MOBY radiometric standards 
are monitored during the year using two smgle<hannel filter radiometers, called Standard Lamp Monitor (SLM) 
radiometers, with filter channels at 412 nm and 870 nm. They are calibrated at NIST on a regular basis. In addition, 
iTansfer*standard artifacts are deployed from NIST to the MOBY field laboratory on a yearly basis to validate the 
calibration protocols and assess the accuracy and stability of the MOBY spectral radiance sources. 

The generally stated uncertainty goal for ocean*color radio mefry is 5 % (ft = 1) for which results in an 
uncertainty of 35 % for chlorophylls? (Hooker e / »?/. 1993). This requirement for high accuracy places stringent 
demands on the uncertainties in water* leaving radiance measurerr>enls with MOBY The uncertainties of MOBY 
mea5urerr>ents of water*leaving radiance L, — without accounting for the effects of stray light — arc estimated to 
be in the range from 4 % to 8 % (ft = I ) (Clark ei al 2002a). 

Spectrogr^hs are imaging systems with dispersive elements and multi*elerr>ent detectors that enable 
simultaneous acquisition of an enure spectrum over some finite spectral width. There are intrinsic limitations in the 
background signal originating from radiation scattered from imperfections in the ^lical elements in the instrument 
This unwanted background radiation, called stray light, while small — on the order of 0 0 1 % or less of the incident 
spectral radiance in a single grating spectrograph — can give rise to unforeseen errors, often much larger than 
anticipated, when the spectral distribution of a source being measured differs significantly from the spectral 
distribution of the calibration source Such a situation is routinely encountered in oceanographic measurements, 
where instruments arc calibrated against incandescent sources with a peak radiance in the shor1*wave infrared 
1000 nm) and subsequently measure the radiance of the ocean, which peaks in the blue spectral region (-*450 nm) 
(Fig. 5.1). Consequently, stray light, if not properly accounted for. can have a significant effect in ocean-color 
research (Clark ei al. 2002a). 

Measurements with the two spectrographs in the MOS systems can be compared in the spectral interval from 
about 580 nm to 630 nm In this spectral range the two spectrographs give different values for or at a comrr>on 
wavelength In Fig. 5.2, is shown for the Top (Top), Middle (Mid) and Bottom (Bot) MOBY arms. The 
difference in the measured radiance in the overlap region is a function of depth, increasing for de^er-lying MOBY 
arms. While the signal is small in this region (af^roximately 1 % of the peak radiance), the lack of agreement in the 
overlap region and particularly its dependence on the spectral distribution of the upwelling radiance are possible 
manifestations of stray light in the system. 

To look for stray light in MOS, the response of the two spectrographs in the MOS205 system to monochromatic 
laser excitation was rr>easured. In this experiment, lasers were directed into an integrating sphere and the radiance 
was measured with each MOBY arm For monochromatic radiation, the entrance slit is spatially imaged on the 
detector. Ideally, no radiation falls on detector elements outside the image In practice, the image is modified by 
scattered light within the spectrograph and every element in the array can in principle have a finite response to this 
rr>onochromalic radiation. Results for the two spectrographs are shown in Fig. S 3. The spectra are similar for both 
spectrogr^hs. There are three components to the image' a strong, sharp peak corresponding to the image of the 
spectrograph entrance slit on the CCD; a broad, peaked structure around the slit image, and a non*zero constant 
component. These three components are similar to specular, haze and diffuse components of reflectance The 
specular component corresponds to the properly imaged radiation; the haze and diffuse components arise from light 
scattered in the spectrograph, principally from the grating. These two components are analogous to the spectral out* 
of'band features commonly ^served in a filler radiorr>eter 

In certain spectral regions, additional features are observed in the CCD response. From physical examination of 
the spectrographs, the additional peaks are associated with a higher order diffraction peak from the grating Optical 
ray-trace modeling of the system confirms this hypothesis Representative MOBY images showing this reflection 
peak in the system response are given in Fig 5 4. Note, for example, the secondary reflection peak centered about 
pixel 1 60 (pixel 725) for light imaged onto pixel 195 (pixel 760) for the BSG (RSG). 

This scattered light causes unknown but potentially significant errors in the measured up-welling radiance. 
Consequently, the MOS instrument was characterized for its stray-light response and the results incorporated into an 
algorithm that was developed to correct its rr>easureTT>enls for stray light In Section S2, we describe the iterative 
^proach used to correct a MOS spectrum for stray light. The derivation of the instrument’s stray-light parameters 
is described in Section S.3 and applied to MOBY data sots in Section S.4. 
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Figuie 5.1 Typical spectral radiances of a lamp'based calibration source and in*water, up*welbng radiance 
at the MOBY site. 



Figuie 5.2. MOBY L„ (typical) for the three MOBY arms showing the lack of agreement in the overlap 
region where both spectrographs make measurements 
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Figure 5.3 MOS205 response lo monoehromauc excl 
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Figure 5.4. MOS205 response (o monochromalic excitation showing reflection peaks. 


5.2 STRAY-LIGHT CORRECTION ALGORITHM 

Each MOS spectrograph uses a 5 12 by 5 12 element CCD array to detect incident radiation When acquiring an 
image, the signals from the central 384 pixels in each column are averaged and a 512-element array is generated 
from each spectrograph. The total signal frwi element / of a MOS CCD is given by the equation 

(5.1} 

where r/A) is the spectral responsivity of element / and L(Af is the spectral radiance of the source being measured. 
Note that r/Aj is the spectral responsivity of element / when considered as part of the spectrograph and includes 
effects such as grating diffraction efficiency and mirror losses 

For monochromatic radiation, an image of the entrance slit is formed on the detector Expressed as a function 
of wavelength rather than array element, this normalized spatial image function is known as the instrument's slit* 
scatter function rT;(X,'-A.) (Kostkowski 1979), with the 'exit slit' in this case determined by the array element's 

spatial width and Aj is the wavelength of element ;'s maximum responsivity. There is a fixed relationship between 
the excitation wavelength and the position of the image on the array As the wavelength changes, the spatial image 
moves across the array. 

Knowing this relationship enables us to determine the fraction of incident light at some wavelength that is 
scattered onto a particular elerrwnt. For example, for the wavelength in Fig 5.3, 0.01 % of the light irrtaged on 
element 420 (BSG) is scattered onto element 300. Assuming each element in the detector array has the same 
average spectral responsivity, the signal from element 300 to radiation at wavelength Xg, is 0.01 % of the signal from 
element 420 

Following Kostkowski (Kostkowski 1979). the total responsivity of element / can be approximated by the 
product of the slit*scaftering function and the maximum responsivity of each array element ^ ) (Fig. 5.5) 
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r;(J.)=^(>.,)a,(X,-X). (5 2) 

Knowing ihe sUi'^atter function, we have a set of 512 coupled equations 

S, =i;(>.,)|cj,(X,->,)i(X)rf>, (5.3) 

that can be directly solved for either or (provided the other variable is known) using a generalized 

least'squares approach (Shumaker 1979) Approximate inverse solutions using Singular*Value*Decomposition 
(SVD) of the applicable matrices are robust and should be directly applicable to this system (Mueller 2002). We 
used a simpler, more intuitive, iterative approach utilizing the inherently discrete nature of the CCD response. 
Separating Eq (5.1) into an in*band {ib) and an out*of*band (oo/i) component and assuming that the source radiance 
IS ^proximately constant over the in*band spectral width, Eq. (S 1) can be written as 

where A, is the wavelength corresponding to the peak responsivity of element /, is the source radiance at ^ , 

and ) IS the integrated in*band responsivity of element / 

15 - 5 ) 

Measuring a calibration source of known spectral radiance, and solving for in Eq. (5.4), the 

integrated in*band responsivity of pixel i can be written 

5 . n{K) 




(5 6} 


The second term on the right hand side of Eq. (S.6) is the stray*light contribution to the total responsivity. 
Utilisng the discrete nature of the spectrograph detector array and substituting Eq (5.5) for ). a recursion 
relation is developed for ^ ) 


) = — 5 




£r;-">(XJ<,,(X-XJ£(XJaX, (5.7) 


where AX is the pixel *to*piKel wavelength spacing and j extends over all 512 elements of the array The original 
input values to the responsivity are the signals divided by the radiance of the calibration source 




s 




(5 8} 


Following the above discussion, the radiance of a source with an unknown spectral distribution can also be 
corrected for scay light. In this case, the relation is given by the expression 

5 , 1 




UK) R,(X,)^a,{X-X)dX • 




where 








(5.10) 


Note that the stray*light corrected responsivity must be used. A stable solution, with changes less than 0. 1 %, is 
reached after 4 to 5 iterations. The solution is robust and insensitive to small variations in the input parameters. 
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Figure 5.5. (a) A spectrograph's absolute spectral responsivlfy, (b) Sllt*scaner function 

for element 300, ^30o('^300 ” Absolute spectral responsivity of element 300. 
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5.3 DERIVATION OF MOBY SLC MODEL PARAMETERS 

There are a lotal of 3 MOS instniments actively used by the MOBY project MOS202 is a 5tand*alone system, 
known aa the MOS Profiler l( )s a ship*board in5tmrr>en( used for bio*optical algorithm devel^meni and validation 
and IS not directly associated with the MOBY buoy MOS204 is used as the sensor in even*buoy MOBY 
deployments; MOS20S is used in odd*buoy deployments. In principle, model shl*scatler parameters need to be 
independently developed for each MOS and for each MOBY deployment, with s^arate factors for each MOBY 
arm, each port (Ea or L„), for the top E, port and for the bottom MOS port. 

Examination of the system response to monochromatic excitation showed that the different arms (both £ and L) 
within an individual buoy had the same slit^carter function. The normalized responses of the MOBY top, middle 
and bottom arms with MOS205 to monochrorrtatic excitation are shown in Figs 5.6 and 5.7 Measurements also 
confirmed that the slit*scatter function remained the same before and after deployment for a particular buoy as well 
as from buoy to buoy for a particular MOS system. However, the slit*scatler function differs between MOS 
instruments (Fig. 5.8). Therefore, we only needed to derive model slit-scatter function parameters for a single port 
on a single MOBY arm for a particular MOS The results were then extended to every MOBY arm, both £, and £j 
ports, and the £, collector for a particular MOS system. Radiation entering the MOS directly from the bottom 
radiance port is imaged differently (Fig. 5.9) onto the CCD array. Consequently, the slit-scatter function parameters 
for the bottom MOS port on MOBY are considered separately and are not addressed in this work. 


Experimental 

In a spectroradiometer equipped with a single element detector, the dispersive element, eg the grating, is 
rotated while the detector is ftxed. The center wavelength of the radiation imaged on the exit slit is a function of the 
angle of rotation of the grating. The instrument’s slit*scatler function can be determined using a single fixed 
wavelength source and changing the wavelength of the radiation centered on the exit slit by rotating the grating in 
small steps. 

In a spectrograph, the grating is fixed and the spectral selection determined by the image of the entrance slit 
onto a reference plane where the multielement sensor is placed as well as the size of the individual elements in the 
sensor array. In this case, if we want to determine the system’s spectral response over some finite spectral range, we 
need a tunable, monochromatic excitation system. NIST recently established a tunable*laser*based radiorr>etric 
calibration facility for Spectral Irradiance and Radiance responsivity Calibrations using Uniform Sources, known as 
SIRCUS (Brown e/ al 2000) In the SIRCUS focility. tunable, monochrorrutic, high-power lasers are fiber coupled 
into an integrating sphere, producing a nearly Lambertian, monochrorruiic source that overfills the entrance pupil of 
imaging systems such as MOS. The radiance of the source is determined at any wavelength by an irradiance 
standard detector, which is in turn directly traceable to primary radiometric standards maintained at NIST. A bench 
prototype of the MOS system was characterized and calibrated in this facility (Habauzit et al 2003). The results 
demonstrated the validity of the stray*light correction approach 

The imaging onto the MOS CCDs in MOBY is slightly different than in the MOS system itself (Fig. 5 9). 
Consequently, a piece *wise character! zaiion approach would not work and the entire MOBY system had to be 
characterized. Because of the logistical problems associated with transporting an object as large as MOBY to NIST 
and to the SIRCUS facility, as well as time constraints imposed by the operational protocols developed for MOBY, 
it was necessary to characterize and calibrate the buoys in situ at the calibration focility at the University of Hawaii 
Marine Center, Honolulu, Hawaii. 

A fiber*coupled, laser*based integrating sphere source (ISS), similar to the radiance source in the SIRCUS 
facility, was developed for use at the MOBY field site The system consisted of a 532 nm Nd Vanadate laser 
pumping a Tiisapphire laser for tunable radiation from 700 nm to 1000 nm or a dye laser with Rhodamine6G dye for 
Tunable radiation from 565 nm to 615 nm The Tiisapphire laser was modified to enable intra*cavify doubling with 
doubled output from 4 1 0 nm to 440 nm. An etalon was placed in the laser cavity to ensure that the output was 
narrow'band (<0.01 nm). Fixed*frequency HeNe, Argon ion, and diode lasers were used for additional 
characterization at discrete wavelengths. 
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Figure 5.7 Normalized MOS205/MOBY2I9 response to monochromatic excitation ftom the three arms 
showing reflection peaks. 
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Figure 5.8. Comparison of MOS204/MOBY2 1 8 (lighJ solid line) and MOS205/MOBY2I9 (dark solid 
line) radiance ports' response (o monochromaiic exciialion 



Figure 5.9. Comparison of MOBY radiance arm (light solid line) and MOS (dark solid line) response to 
monochromatic excitation. 
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Figure 5.10; Scbcmniic of cxpeunental setup at the MOBY field site ai the Univeisity of Hawaii Marine 
Center, Honolulu, Hawaii, (a) Exterior setup showmg location of lasers and calibration tcni. (b) 
experimental setup in the calibration lent for mcasunng the ISS radittncc with a coUxior on the middle 
MOBY arm. 
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The output from ihe lasers was sent into a laser power*controller/siabllizer to reduce power fluctuations in the 
beam. Portions of the beam were also sent into a wavemeter to determine the wavelength of the laser radiation and 
into a Fabry*Perot spectrum analyzer so that the mode structure of the laser output could be monitored. The lasers 
were kept in an air*conditioned sea van and radiation was f)ber*coupled into an mtegrahng sphere in the MOBY 
tent A portion of the optical fiber was placed in an ultrasonic bath to reduce effects of speckle in the measurement 
of the output of the ISS. From the sea van, the optical fiber was placed in a garden hose for protection and through 
an electrical conduit to the MOBY tent The experimental setup is shown in Fig. 5.10. 

The sphere radiance was determined using a calibrated irradiance Dansfer*.'<tandard radiometer. The radiometer 
is a six*element silicon*transmission trap detector (Eppeldauer and Lynch 2000) It was calibrated for spectral* 
power responsivity by comparison against a cryogenic radiometer and was equipped with a precision aperture of 
known area. The sphere radiance is given by the ratio of the measured signal to the trap responsivity times a 
geometrical flux transfer constant G|> 


with <Jo given by 



_ 


(5.11) 



(5.12) 


where we have +d'^ , Ai (r,) is the sphere exit port area (radius), Ai (r^ is the trap aperture area 

(radius), and d is the distance between the two apertures 


The sphere was equif^ed with a silicon monitor photodiode During MOBY measurements, the monitor was 
calibrated for the sphere radiance over a spectral range of interest using the reference standard trap detector, and 
then MOBY measurements were made The monitor signal was recorded for each MOBY measurement. 


Communication was established between the laser sea van and the MOBY tent using wa)kie*talkies or a marine 
band radio. It was important to properly couple the laser radiation into the optical fiber. During this alignment 
process, the monitor signal was continuously displayed on a computer in the MOBY tent. Coupling of the laser light 
into the fiber was adjusted, and the effect on the monitor signal was noted The laser was properly coupled into the 
fiber when the monitor signal was maximized. The monitor signal was dLsplayed on a computer screen in the 
calibration tent. Using commercially available software, the screen was also displayed on a computer in the sea van 
to simplify the fiber alignment process. 


Sht'seatter parameters were determined for the BSG and the RSG coupled to the MOBY arms. To properly 
model the slit*scaftei function (SSF) for the stray* light correction (SLC) algorithm we need the following 
parameters (a) the single pixel in*band area; (b) a fit to the OOB slit*scatier function; (c) model for the reflecbon 
peaks and (d) model for addihonal light scatted onto the BSG CCD from radiation in the spectral range ftom 640 nm 
to 690 nm In the following subsections, the basic steps undertaken to derive mode) parameters for the slit*scatter 
function are outlined for MOS205/MOBY 219 middle arm L„ port for the BSG and RSG, respechvely. 


MOS205/MOBY2I9 BSG fn-band area 

The stiay*hght correction terms on the nght hand side of Eqs (5.7) and (5.9) are divided by the normalized in* 
band area, N. N is given by Eq. 3.3, where the integration is over a 'suitably defined' in*band area 

= (5.13) 

To determine the in*band area, the laser output was changed in small increments (0. 1 nm to 0.2 nm) from 590 nm to 
600 nm and an image was acquired at each wavelength. The raw data (ADU/pixel/s) were divided by the sphere 
radiance, then normalized to a pixel's maximum responsivity. The wavelength of each pixel's maximum 
responsivity was subtracted, leaving each pixel's normalized responsivity, centered about a relahve wavelength of 0 
nm, for three adjacent pixels, P420*F422. Fig. 5.11 shows the normalized in*band responsivity of these pixels. 
Numerically integrating the line shapes gives in*band areas for ihe three pixels that agree to within 1%. Changing 
the limits of integration from ± 4 nm to ± 6 nm changed the areas by less than 0.2 %. Fitting to an analytical 
function also gave results in agreement with the numerical integration The imband area and associated relative 
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siandard uncertainty are listed in Table 5. 1 . ln*band scans were repeated over the range from 420 nm to 430 nm to 
delennme the spectral dependence of the bandwidth No dependence was observed 


Table 5.1 . Values and uncertainties in the sliKcarter function derived uncertainties for 
MOS205/MOBY2 1 9 Lu middle arm. 


Parameter 

Variable 

MOS2 05/MOBY 
BSG 

Unccrtaintv 

\%] 

MOS205/MOBY 

RSG 

Uncertainty 

I%1 

Normalization 

•V 

0 9fi5 

1 

1 137 

1 

Bandwidth 

bw 

0.5S05 

0.25 

0.8047 

0.25 

OOB SSF 

m, 

0.1142 

2 

0.164 

2 



1 13.17 

2 

9042 

2 


m, 

I.l25e*5 

25 

1.5e*5 

25 

Reflection Peak 
Position 


92 4S5 

0.25 

343 0125 

0.25 


h 

0.6m 

0.25 

0.5859 

0.25 

Reflection Peak 
Area (P<I00} 

Ot 

0.02599 

3 

0.0279 

3 


Q; 

295 0 

3 

713 997 

3 



33.0 

3 

25.4843 

3 

Area I P> 100} 

h, 

0.02599 

3 




ft. 

346.2642 

3 




hi 

13.2884 

3 



Off CCD 
Amplitude 

ScaleFactor 

0 14 

10 

NOP 


Reflection Peak 
Dispersion 

Ct 

0.6108 

0.25 

0.5859 

0.25 


MOS205/MOB ¥2i 9 BSO out-of-hand slii-scarter futtciion component 

In modeling the slit*scatter function, we are considering the normalized relative response of the array to 
rr>onochromalic excitation. To determine the out*of*band (GOB) component of the slit* scatter function, the laser 
excitation wavelength was fixed at the peak of the response of one of the in*band pixels and a number of images 
were acquired and averaged. The OOB slit*scarter function is shown in Fig. 5. 12. The spectrum, with the center ± 5 
pixels masked, was hi to a double exponential of the form 

V = m, exp(-2|jr-x,|)/m,)+m, , (5.14} 

where x = column number and is the column number with the peak response to the excitation wavelength, and the 
m,’s are constants determined from the fit The fit is shown by the solid line in Fig. 5.12, and hi parameters are 
listed in Table 5. 1 . Note that there is an apparent 2*pixel shih in the center pixel between the imaged peak (centered 
at column 419} and the fit peak =421). This shih is ignored in the model of the slit* scatter function 


MOS205/MOBY219 BSG Reflection peak;: 

The amount of light scattered onto the CCD in this reflection peak is significant (I is therefore necessary to 
include a model for these peaks in the general model for the slit*scatter function For the model, we want to know 
what fraction of the incident light imaged onto some pixel is measured by the pixel under consideration. The 
rr>odeling is complicated by the fact that the reflection peak area is a strong function of position and the rale at which 
It moves across the CCD with changing wavelength differs from the primary image. 

To model the incident light scattered onto a particular pixel, we first need the relationship between the reflection 
peak position and the excitation peak position. This relationship, shown in Fig 5 13, was fit to a linear regression, 
given by the solid line in the figure. Values are given in Table 5.1 . 

To calculate the fraction of the incident light that is in the reflection peak, we want to compare the integrated 
areas of the reflection peaks to the areas of the fundamental peak. The primary peaks in the data set are first 
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normalized (o 1. The relative reflection peak data, plotted a.< a function of separation from the primary image (in 
pixels), are shown in Fig 5 14. There is a large variance in both the amplitude and the total area under each 
reflection peak, but we need to consider that not all excitation wavelengths were centered on a pixel’s response, so 
the reflection peaks may not be properly normalized to a pixel’s peak response. Consequently, it was nece?isaiy to 
correct the amplitude of the reflection peaks in Fig. 5 14 for the excitation offset. 

The in*band data set was fit to a Gaussian distribution given by 


(5.15) 

where the p,*s are constants determined from the fit. 

We get a reasonable fit, shown in Fig. 5.15. though the area is slightly ur>derestimated because the fit does not 
represent the righl*hand shoulder well. The calculated area from the Gaussian fit agrees with the numerically 
integrated value to within 1 %. In Eq. (S.15), pi is the amplitude; p; is the center or offset, and pj is the width To 
correct the excitation peak, we fit the primary image for the reflection peak data set to a Gaussian lineshape with a 
fixed width determined by the in-band fit Note that to obtain the proper width for the fitting routine, it is necessary 
to convert the width from wavelength to pixel by dividing the width (nm) determined by the fit in Fig 5.15 by the 
pixel'tO'pixel spacing of the BSG (0.58 nm/pixel). 

With the width fixed, we are only concerned with the amplitude determination (or the value of pi) in Eq. 3.4. 
Each primary image peak was fit to Eq. (5.1 S), and results are summarized in Table 5.1. Corrections as large as 
26 % are applied to the reflectance peak amplitude data set. Note that there is o^en a small shift in the peak position 
as well The shift is less than O.S pixels and is not included in the modeled reflection peak position. 

The corrected, normalized reflection peak data set is shown in Fig. 5. 1 6. The 440 nm data set had an additional 
correction applied to make the haze component similar to the other data sets. One possible explanation for the 
^parent increase in scattered light is a wavelength*instability problem with the 440 nm diode laser used to acquire 
the specTTum This would broaden the primary image and lead to an overestimation of the reflection peak amplitude 
in the modeled slit*scatter function 





For the .vlit*scatier function model, we are interested in the fraction of the in*band radiation .scattered into the 
reflection peak. The amplitude*corrected reflection peaks are fit to Gaussian line shapes, giving the reflection peak 
area as a function of excitation wavelength or peak pixel. The reflection peak area as a function of position is shown 
in Fig. 5 17. The solid line is the fit to the data using Eq. (5.16) for pixels greater than 1 00 and Eq. (5. 17) for pixels 
less than 100, ie. 


y=<j,/()+(x/o.)“‘).and (5.16) 

V = /(I +((350 -z)/^ )'’•)• (5.17) 

where x refers to the pixel under consideration. 

The area must be further corrected for the increased rate at which the reflection peak moves through a particular 
pixel compared with the primary image. The primary image moves at a rate of approximately 0 58 nm/pixel while 
the reflection peak's dispersion is 0.3S nrrVpixel. This will be discussed further, in subsections to follow, with 
regard to the red spectrograph. 
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RerectMn Peak PesAlon (Pfcrol) 

Figure 5.13 Location of ihe MOS205 BSG primary image versus the position of the reflection peak, 
with the solid line a fit ro the data 



Figure 5.14 MOS205/WOBY2I9 L„Mid BSG reflection peak data set, uncorrected for the amplitude 
of the excitation peak 
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Figure 5.15 Gaussian fit to the MOS205/MOBY219 BSG in*band data set. 



Figure 5.16' MOS205/MOBY2I9 LuMid BSG reflection peak data set corrected for the 
amplitude of the excitation peak 
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Rerectksn Peak Poslllon (Plxal) 

Figure 5.) 7' Reflection peak area as a function of position; tPe solid line Is a fii to the dau. 
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Figure 5.18; Response of BSG pixels 400 lo 51 1 and full RSG response for discrete laser excitation in 
the range from 640 nm to 690 nm 

M205/MOBY2I9 BSG Of -CCD Scartenng 

There is additional scattering onto the right edge of the CCD from light not directly imaged onto it In Fig. 
5. 1 8, we show the responsivity of pixels 400 to 5 1 1 for the BSG from discrete laser excitation in the range from 640 
nm to 690 nm. Tunable lasers in the 640 nm to 690 nm spectral range were not available during the Reid 
mea5urerr>ents and we do not have enough spectral data to directly model BSG stray light effects in this range 
Consequently, the magnitude of this contribution was empirically fit to validation source data sets (discussed in 
Section 5) by including a quadratic term in the BSG slii*scaiter function in the overlap region (pixels 400 to 511). 

MOS205/MOBY2J9 Red Spectrograph (RSG) 

The same set of model sUi^catter parameters was determined for the MOS205/MOBY2 1 9 RSG, with the 
exception of the off CCD term. Flux transmitted through the dichroic beamsplitter, incident on the RSG, falls off 
rapidly below 600 nm. Essentially no light below 550 nm is incident on the CCD, so a term to account for off*CCD 
scattering is not included in the RSG model slibscatler function. 


MOS20S/M2I9RSG in-band area 

Images were acquired as the excitation wavelength was tuned in small steps over the spectral range from 740 
nm to 750 nm The in*band areas of three adjacent pixels P756, P757, and P758, shown in Fig. 5.19, were 
determined. The areas agreed to within 0.25 %; changing the limits of integration from ± 4 nm to ± 6 nm changed 
the in*band area less than OJ %. The in*band area and uncertainty are listed in Table 5. 1 . 


MOS205/M2I9 LuMid RSG OOB Slii-sca/ter Fuftction 
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The peak •normalized GOB slii-scaiter function, along with the double^exponeniial fit to the data (Eq. (5.13)] 
aie shown in Fig 5.20. Nore that the reflection peak was marked for this fit as well as the in*band ± 5 center pixels. 
There is a 06 pixel offset from center between the image and the fit. This small offset is ignored in the model. 
Parameters are listed in Table 5 1. 

MOS205/M2J9LuMid RSCs ReJIecrion peakf 

Images were acquired for excitation wavelengths ranging ^m 640 nm to !^00 nm. A linear fit to the 
relationship between the reflection peak and the excitation peak positions is shown by the solid line in Fig 5.21 
Next the central peaks were fit to a Gaussian lineshape, with a width determined from the in*band scan. As with the 
BSG, these data are used to properly normalize the reflection peak amplitude For the RSG, the dispersion is 0.SO4 
nnv'pixel Uncorrected and corrected reflection peak data sets are shown in Figs. 5.22 and 5.23. Note that there is 
an apparent increase in the amplitude of the reflection peak for the first 3 wavelengths, where diode lasers were used 
for excitation As with the BSG 440 nm data set, it may be a consequence of the multi*mode spectral output of the 
diode lasers used for this spectral region. 

The amplirude<orrected reflection peaks are fit to Gaussian functions with a linear offset, giving the reflection 
peak area as a function of excitation peak pixel. The areas are again corrected for the rate at which the reflection 
peak rr>oves through a particular pixel. For example, consider the RSG fine scans arour>d 740 nm. In this spectral 
region, a reflectance peak is observed around 700 nm (pixel 702). as shown in Fig. 5.24. As the excitation is 
scanned in srrtall steps, the response of pixel 702 can be measured and compared with the spectral width of the 
reflection peak in Fig. 5.24 Fig. 5 25(a) shows the reflection peak observed for 740.3 nm excitation, while Fig. 
5.25(b) shows the spectral response of pixel P702. Both are plotted in wavelength scale using a dispersion of 0.SO4 
nm/pixel for Fig. 5.25(a). Note the reduction in the width of the P702 reflection peak with respect to the reflection 
peak image from monochromatic excitation. This apparent reduction occurs because of the difference in the spectral 
dispersion of the primary and reflection peak images, and must be taken into account in the rr>odel of the reflection 
peak amplitude 

The final results are graphically illustrated in Fig 5.26 with the fit to Eq. (5.16) shown by the solid line. Fit 
parameters are listed in Table 5.1 
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Figure 5.19: CeniCT*wavelength sublrscted, peak normalized. MOS205/MOBY219 RSG in*band 
scans for pixels 756.757, and 758. 



Relallvo Column (Pixel) 


Figure 5.20 Peak normalized MOS205/MOBY2 1 9 RSG ou(*of»band slit*seatter function. 



Figure 5 J I . MOS205/MOBY2 19 oxcitaiion peak position as a function of reflection peak position. 
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Figure 5.23: MOS205/MOBY2 19 corrected reflection peak data set. 
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Figure 52A: Relative RSO response to monochromatic excitation at 740.3 nm. 



Wavelength (nm) 

Figure 5.25: (a) Reflection peak for RSG response to monochromatic excitation at 740.3 nm 
(b) Response of pixel 702 as the excitation wavelength is changed from 695 nm to 705 nm. 


109 







Ocean Opuc< Proiocols For Saielliie Ocean Color Sensor Validaiiorw Re\‘ision 4, Volume VI 



Figure 5.26: MOS205/MOBY219 reflection peak amplitude as a function of reflection peak 
position 


5.4 STRAY-LIGHT CORRECTION OF MOBY 

The parameiers developed in Section 5.3 for MOS205/MOBY2I9 were incorporaied into ihe scay-light 
correction algorithm described in Section 5.2. The algorithm was used to correct MOBY calibration and in-water 
data sets MOBY was calibrated using a lamp*illuminated integrating sphere calibrated for spectral radiance. The 
iterative stray-light correction procedure described in Section S 2 [Eq (5.7)] was applied to the measured signal. 
Iterations were repeated until a stable solution was reached; a solution that changed by less than 0.1 % from the 
previous iteration was considered stable. Typically, a stable solution was reached after 4 to 5 iterations. The 
uncorrected responsivity and the stray*lighr corrected responsivities are shown in Fig. 5 27. In Fig. 5J8, the ratio of 
the corrected responsivity to the uncorrected responsivity is shown. Stray-light is responsible for approximately 
1 0 % of the signal over much of the spectral range, increasing rapidly below 400 nm for the BSG. and below 
640 nm for the RSG. The rectangular increase in the corrected over uncorrected ratio (reduction in stray-light) for 
the BSG near 500 nm occurs in the region where the reflection peak coincides with the primary image 

Using the stray-light*corrected responsivities, MOBY up-welling radiance data are corrected using Eq. (5 9). 
Uncorrected and corrected radiance measured by each arm is shown in Fig. 5.29, and the ratio of uncorrected to 
corrected up-welling radiance in Fig 5.30. Note the dramatic improvement in the agreement between the two 
spectrographs in their overlap region in the stray-light-corrected results. There is also a significant increase in the 
measured up-welling radiance below 400 nm 

MOBY is used to vicariously calibrate SeaWiFS by providing measurements of water-leaving radiance. The 
water-leaving radiance is determined ^m the up-welling radiance according to Clark (Clark ei al. 2002a, Clark ei 
al. 2002b). The water-leaving radiance was calculated for stray-light-corrected and uncorrected data sets The 
correction factor — or ratio of corrected to uncorrected water-leaving radiance — is shown in Fig. 5.3 1 . We again 
see a strong increase in the corrected water-leaving radiance below 400 nm, reaching a value of ^proximately 40 % 
for 360 nm radiation. The corrected water-leaving radiance is less than the uncorrected values for wavelengths 
longer than 500 nm The decrease rises to approximately 20 % for the BSG beyond 600 nm 
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Band* integrated siray*light*corTection factors for water*leaving radiance for the SoaWiFS bands are listed in 
Table 5.2. The correction Actors range from +8.5 % for Band I to -2.5 % for Band 4 


Table 5.2. Band*integratcd stray*light*coireclion factors for 
water*leaving radiance in the SeaWiFS bands. 


SeaWiFS 

Ratio 

Corrected^ncorrected 

Band 

L. 

1 

1.085 

2 

1.033 

3 

1.012 

4 

0.976 

5 

0.979 

6 

1.014 
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Waveiengih (nm) 

Figure 5.27' Uncorrecled (diamonds) and stray^lighUcorrecied (solid line) responsivity of 
MOS205/MOBY2I9 radiance arms 



Wa%elen^li (ntn) 

Figure 5.28: Ralio of Ihe MOS205/MOBY2I9 UMid radiance arm eorrecled responsivity lo the 

uncorrecled responsivity. 
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Figure 5,29: Uncorrecled (thin solid line) and corrected (thick solid line) up-welling 
radiance measured by each MOS205/MOBY219 arm, (Bottom panel) Expanded view 
in the UV region. 
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Figure 5.30; Ratio of corrected to uncorrecled up*wellmg radiance measured by each 
MOS205/MOBY2I9 arm. 



Figure 5.31 Ratio of 8lray*lighi*corTecied to uncorrecled water*leaving radiance measured by 
MOS205/MOBY2I9 
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5.5 ALGORITHM VALIDATION 

To validate the 8tray*light correction mudol, a Colored Source (CS) was developed. Different relative specoal 
distributions were realized using glass and loterference filters and a lamp'illuminatcd integrating sphere source The 
spectral radiance for the CS configurations was determined separately using a calibrated specoo radiometer The 
SLC algorithm was validated by comparing the corrected spectral radiances from measuicmcots with MOBY to (be 
blown CS spectral radiances. 

A lamp'illuminaled integrating sphere source, model OL420 from Oplionic Laboratories,'*' was used for the CS. 
Briefly, it is a 20.32 cm diameter, externally iDumJnaied sphere that is coated with barium sulfate. The exit aperture 
IS 5 72 cm in diameter and the entrance aperture is orientated 90'’ with respect to the exit aperture. Quartz windows 
on both apertures protect the sphere from contamination. A monitor dclecior. fitted with a pholopic filler, is 
mouoicd in the wall of the sphere and its output voltage is recorded during operatioru 



Wavelrogib [nm) 


Figure 5.32: Spectral radiance of the OL420 Colored Source with different filters. The 
radiance of the CS with the 1400 nm interference filler, which is a very low level, is shown on 
the righl'hand axis VXR measurements of the CS are given by the open symbols 


For Ihe CS, we utilized the filter slide that is located between the lamp and the eutranee aperture. Six CS 
configurations were possible, with each filter in a corresponding filler slide for ease of use and minimization of 
con lamination through handling. Three colored glass filters were used: Schofl BG28 and BG39^', and CVf'* PER 
filter, producing blue, turquoise and photopic (grecc) spectral distributions Three ioterfenTnee filters were used, 
centered ai 400 nm, 700 nm, and 750 nm, with full* width half* maximum (FWHM) bandwidths of 10 run, 45 nm and 
1 0 nm, respectively. 


Certain commercial oqurpmenl, instruments, or materials are ideniified in this paper to foster understanding. Such 
identification does not imply recommendotiou or endorsement by the National (nsntute of Standards and 
Technology, nor docs it imply that the materials or equipment identified are necessarily the best available for the 
purpose 


i; 


Schott Glass Corporation (www.us.sehoft.com). 
CVf Laser Corporation. Albuquerque, NM. 
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Figure 5.33: Uneoirecled (Ihin line) and stray^lighl-corrected (dark line) and ‘true’ (gray circles) 
measured radiance of the Colored Source with (a) BG28, (b) PER, and (e) [750 fillers. 
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The OL420 was calibrated for spectral radiance on the NIST Facility for Spectroradiometnc Calibrations 
(PASCAL) (Walker ei al. 1987) in March 2001. The .spectral radiances for the PER, B028, and BG39 
configurations of the CS were determined at NIST m May 2001 using a reference spectroradiometer The 
instrument is a double*grating, scanning monochromator (Brown and Ohno 1999). It has been fully characterized, 
with a stray* light factor of less than 10*". For the CS calibration, the bandpass was set to 5 nm FWHM and a 
wavelength step of 5 nm was used. For determinations of the CS ^ectral radiances, the unfiltered OL420 was used 
to calibrate the scanning monochromator. Typical spectral radiances for the six configurations are shown in 
Fig. 5 32. The reference spectroradiometer is not a field instrument, and was not transported to the MOBY field 
site. 

During the field deployments, the EOS/NfST Visible Transfer Radiometer (VXR) was used to verify the CS 
spectral radiances. The VXR is a siK*channel filter radiometer (Butler er al 1 999), and comparison to the reference 
spectroradiometer results could only be done where the VXR filter channels coincided with measurable output for 
the CS configurations. The agreement with the PER was *0.8 % at 548 nm and +1 9 % at 6<> I nm; with the BG28, it 
was *3 % at 44 1 nm and *0.7 % at 54K nm, where negative values indicate the VXR values were lower than the 
reference spectroradiometer. The VXR results are shown as symbols in Fig. 5.32. 

The CS was measured by the radiance pons on the upper or middle buoy arms. The results for the source 
BG28, PER, and [750 filters are shown in Fig 5.33. The spectral radiances derived using MOBY are shown wiih 
and without correction for stray light along with the spectra) radiances for the CS configurations determined using 
the reference spectroradiometer. Within the in*band spectral regions, the corrected MOBY values agree with the CS 
spectral radiances to within 5 %. 

The out*of*band regions are an important test of the stray*light*correction algorithm. The BG28 configuration 
decreases by an order of magnitude in the crossover region between the BSG and the RSG, similar to the behavior 
seen in MOBY waters [see Fig. 5.33(a)]. For this configuration, the correction algorithm was sensitive to the 
parameters for the off*CCD scattering, resulting in under* or over*correction in the overlap region We determined 
the off'CCD parameters by minimizing the differences between the two spectrographs for the BG28 configuration 

Figure 5.33(a) shows that MOBY overestimates the spectral radiance as the wavelength decreases if not 
corrected for effects arising from stray light. The PER configuration, shown in Fig 5.33(b), is an excellent test of 
the ability of the SLC algorithm for the shorter measurement wavelengths in the BSG. The spectral radiance in the 
region below 450 nm for the PER configuration decreases by more than IC^ compared to the maximum value 
Uncorrected, the MOBY values are discrepant by more than three decades in this spectral region. With correction 
for stray light, the agreement is in^roved substantially 

The inability to correct at the shortest wavelengths, below 400 nm, is understood in terms of the low system 
responsivity in this spectral region 

A test of the algorithm for the RSG is shown in Fig 5.33(c) using the 1750 interference filter (peak 
transmittance at 750 nm). In the out*of*band region, the uncorrected MOBY values are discrepant by a factor of 10 
to 100, but the corrected values are in good agreement. 


5.6 UNCERTAINTIES 

A Monte Carlo approach was used to estimate the uncertainties in the stray*light*correction Actors for MOBY 
up* welling and water*leaving radiance measurements. Model slit*scafter function parameters were input into the 
stray*light correction algorithm using mean values given in Table 5 1 and a random normal distribution with 
uncertainties given in Table 5.3. Uncertainties are approximately a ^ctor of2 larger than those given in Table 5.1 to 
account for possible additional uncertainties in the model parameters arising fi^m changing environmental 
conditions and other unknown factors. The full sfray*light correction algorithm * correcting both the responsivity 
and the in* water data * was run 100 times for a typical in* water spectrum and mean correction factors and 
uncertainties calculated Typical results for the BSG are shown in Fig. 5.34 for MOS205/MOBY2I9 LuTop Both 
the correction factor and the associated relative uncertainty is less than 1 % over the spectral range from 420 nm to 
590 nm and both increase rapidly below 400 nm and above 600 nm, rising to approxirrtaiely 8 % at the ends of the 
array: 360 nm and 640 nm, respectively. 
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Similarly, the uncertainly In the stray-light'correction factor for derived waler*leaving radiance was delennmed, 
as shown in Fig. S 35. The uncertainty In Lw is larger than the uncertainty in Lu, with a 14 % uncertainty in the 
5Uay*lighl*coiTeclion factor at 360 nm. In Table 5.4, the relative uncertainties In the band*averaged mean stray* 
hghl*eo!Teellon feelers are given for SeaWiFS bands. 


Table 5.3. Uncertainties in sllt*scafter function 
model parameters used in Monte Carlo simulations 
of MOS205 rrwasuremenls. 


Table 5.4: Relative combined standard uncertainties 
in the band*averaged stray*light*correction factor 
for SeaWiFS bands. 


Parameter 

8SG 

Uncertainty [%] 
RSG 

Uncertainty [%] 

Normalization 

2 

2 

OOB SSF 


5 

3 


5 

3 


25 

25 

Reflection Peak Amplitude P<100 


5 

5 


SeaWiFS 

Band 

Ratio CorrectedTincorrected 
Uncertainty 
[%] 


I.0S5 

0.83 

2 

1.033 

0.36 

3 

1.012 

0.15 

4 

0.976 

0.21 

5 

0.979 

0.58 

6 

1.014 

NOP 
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02 

5 

5 

oj 

5 

5 

Amplitude (P>I00) 


a> 

5 

5 

oj 

5 

5 

Off CCD Amplirude 
10 
NOP 


5.7 DISCUSSION 

The SeaWiFS band*averaged waier*leavmg radiance stray*lighl*coiTeciion factors measured by MOBY during a 
rwo-month deployment (M220} is shown In Fig. 5 36. [n Figs. 5.37 and 5.38. the correcied/uncorrecied SeaWIFS 
band*averaged ratios are shown for 4 years of deployments Each grouping in the figure is a separate deployment. 
The correction factors have been stable over the entire deployment sequence, implying that the MOBY imaging and 
the MOS slit'scafter functions, along with the ocean*color, have rerrtalned stable over this tlrr>e frame. The only 
exception was for the first three deployments of the MOS204 sensor (even buoys. Fig. 5.38}, when the Band I 
correcnon fector was 2 % higher than for the other deployments. For these three deployments, low hydroxyl (OH) 
optical fiber was used to couple the radiance heads to the Instrument bay. This fiber transmitted less light than the 
other fibers In the blue spectral region, in particular below 400 nm For these three deployments, then, the relative 
stray*light contribution to the total signal in the blue spectral region was larger than for subsequent buoys, resulting 
in the larger Band I stray*light*correction factor for these buoys The observed radiometric stability of the MOBY 
systems enables us to correct previous deployments for stray light with confidence 
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Mg. 5 37 Band*avcragcd Lw %iray*lijjhl corrcciions for ScaWiFS baud' mcasuicd bv 
odd MOBY buoys from 1998 through 2002 
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Fig. 5 36 Band*averaged Lw stray*)ighi eoTreellon factors lor ScaV^'iFS bands 
measured by MOBY during deployment 220. 
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Impact on Sea WiFS 

MOBY is used for the vicarious calibratioo of SoaWiFS bauds I through 6 (Bamcs et al. 2000). Any changes 
in MOBY'dcnvcd walcr*lcaving radiances are dircclJy reflected in SeaWiFS calibration coefficients for these 
channeU. The effect of the MOBY stray*light correction is lo increase the v*ater*leav'ing radiance of Sea WiFS 
Bonds I through 4, and decrease the radiance of Bands 5 and 6. ScaWiFS uses the Band 2 (443 nm) to Band S (5S5 
run) ratio to radiometncally determine chlorophyll a concentrations in oligotrophic waters (chlorophyll 
concentrations about 0.3 pg^l); the Band 3 (490 nm) to Band S rano for mcsotrophic ualcrs (chlorophyll 
concentrations ranging from 0.3 pg/l to 1,5 pg/l); and the Band 4 (510 nm) to Band 5 ratio for cutrophic waters 
(chlorophyll concentrations > 1.5 pg/'l) (O'Reilly et at. 1998), Corrccring for stray light increases the band ratio in 
oligotrophic waters by 5 % to 6 % and approKiroalcly 3.5 % for waters with higher chlorophyll a concentrations. 
The change in band ratios corresponds to a 25 % lo 35 % decrease in chlorophyll in oligotrophic waters and a IS Va 
decrease in mesotrophic and eutrophic waters. The stray*light correction of MOBY reduces global UKan 
chlorophyll a concentrations measured by SeaWiFS by 15 % to 20 %. When combmed with additional ptoceduraJ 
changes implemented in the 4“ reprocessing of the SeaWiFS data set in July 2002 (SeaWiFS 2002), the result was a 
mean monthly decrease in global chlorophyll a concentranou of 6 %. In addition, the stray* Lght correction greatly 
reduced the magnitude of an issue with SeaWiFS measuring negative waier*lcaving radiance in coastal regions with 
complex aerosol distributions. 
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Figure 5.38: Band*averaged L« stray*light*corTectioQ factors for SeaWiFS bands measured by even 
MOBY buoys from 1 998 through 2002. 

Impticationx for oigorithm devetopmenl 

The derivation of remotc*sensing*based ocean color data products such os chlorophyll a involve the integration 
of measuremcols by (1) the satellite sensor {e.g. SeaWiFS), |2) the vicarious calibratioo sensors (eg MOBY), and 
(3) the instruments used to develop the bio<ipncal algorithms relating the physical propertic« of the ocean (e.g. 
chlorophyll a concentration) to a radiometric measurement Errors io measurements of any one of the three 
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components of the measurement chain will significantly impact the validity of the final data product Instruments 
involved in ocean color measurements are typically calibrated against incandescent sources such as lamp* 
illuminated integrating spheres or reflectance plaques. These types of calibration sources have very different 
spectral distributions than the distribution of uj> welling or waier*leaving radiance and effects of siray light on the 
instrument response need to bo considered for these insfrumonis 

While considerable time and effort can be (and has been) expended to fully characterize and correct a satellire 
sensor or primary vicarious calibration station for stray light, it is much more difficult to extend the detailed 
characterization to the mynad instruments used by different groups to develop bio*optical algorithms. No general 
protocol has been developed to characterize and correct instrumentation used to develop bi^optical algorithms for 
stray light. However, effects of stray light are effeciively canceled when the calibration source distribution 
approximates a measured ocean color spectrum. We have developed a radiometncally stable, spectrally tunable, 
solid'Siate calibration source for ocean color. This unique source approximates the spectral distributions of waters 
with varying chlorophyll concentrations, and may be used to calibrate, characterize, and to correct instruments’ 
responsivibes for stray light. 

We have developed the prototype tunable solid*state source using Light Emitting Diodes (LEDs) (Brown ei al. 
2002) The source has 40 different, independently controlled LED channels comprised of 3 to 6 LEDs (depending 
on their output power) organized into 4 heads that mount onto 5 cm ports in the wall of an integrating sphere The 
flux frw> each channel can be independently changed to model the spectral distribution ofwatei*leaving radiance in 
different chlorophyll concentration waters. The output is monitored using a reference spectroradiometer integrated 
into the system. The source is under computer control and is designed to be autorruiically adjusted lo match any 
desired input target spectrum, within the bandwidth limitations of the individual LEDs (*- 20 nm), to within 
^proximately 5 % over the range from 390 nm to 550 nm 

The source has been tuned to match the water*leaving radiance from blue (oligotrophic), blue*green 
(mesotrophic), and green (eutrophic) waters, covering a wide range of chlor^hyll concentrations The target 
spectral distributions and the solid*staie source output are shown in Figure 5.39. There is reasonable agreement over 
rr>ost of the spectral range The sphere radiance is low below 390 nm. This is partially due to the low refleciance of 
the integrating sphere coating With a new coating, we expect an order of magnitude increase in the radiance below 
400 nm. The source is currently missing diodes al 360 nm, 435 nm, and 455 nm. The radiance around 550 nm is 
low due 10 weak diode emission; a problem thai is easily addressed with the addition of extra 550 nm LEDs. The 
source stability and its applicability to the development of stray* light correction algorithms are currently under 
investigation. 

We propose to correci instrumonis for stray* light errors in water-leaving radiance measurerr>ents using a 
variation on the basic approach developed for SeaWiFS by Wang, ei at. (2001) Knowing the total spectral 
responsivity of the SeaWiFS bands, a stray light, or spectral our*of*band correction, is currently applied to each 
channel based on a semi-analytical water* leaving radiance rr>odel devel^od by Gordon for water with varying 
chlorophyll a concentraiions (Gordon 1995; Wang ct al 2001). In this case, we don’t know the total spectral 
responsivity of the instruments. Instead, we used the solid*stale source to mimic the varying spectral water*leaving 
radiance spectral dLstributions and derive correction factors based on the ratio between the known radiance of the 
source (as measured by our reference spectroradiometer) and the radiance measured by the ocean color in8trurr>enl8. 


5.8 SUMMARY 

A simple iterative procedure was devel^ed to correci MOBY measurements for stray light. Several fully 
operational MOBY systems were characterized in sUu using a tunable laser-based, fiber*coupled iniegiating sphere 
source that was developed and installed at the MOBY field site at University of Hawaii Marine Center, Honolulu, 
Hawaii for this purpose MOBY buoys with MOS204 and MOS205 detector systems were characterized and model 
slit*8catter functions developed. MOBY measurements of a Colored Source validated the stray*light*corrected 
measurerr>ent8. In-water MOBY measurerr>ent8 of and L. were corrected for stray light, and the effect on 
SeaWiFS bands and derived chlorophyll*a concentrations wore evaluated. 
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Figure 5,39: Targol and LED source spcctn for (a) blue, (b) bluc*grecn and (c) green 

waters. 
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